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Abstract 

Spin-density-functional theory (SDFT) is the most widely implemented 
and applied formulation of density-functional theory. However, it is still 
^ ' finding novel applications, and occasionally encounters unexpected prob- 

lems. In this paper we first briefly describe a few of the latter, related to 
\^ ■ issues such as nonuniqueness, noncollinearity, and currents. In the main 

part we then turn to an example of the former, namely SDFT for the 
l/"} , Heisenberg model. It is shown that time-honored concepts of Coulomb 

DFT, such as the local-density approximation, can be applied to this (and 
other) model Hamiltonians, too, once the concept of 'density' has been 
suitably reinterpreted. Local-density-type approximations for the inho- 
mogeneous Heisenberg model are constructed. Numerical applications to 
finite-size and impurity systems demonstrate that DFT is a computation- 
ally efficient and reasonably accurate alternative to conventional methods 
of statistical mechanics for the Heisenberg model. 
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1 Introduction 



Density-functional theory (DFT) is widely recognized as a versatile and pow- 
erful tool for practical calculations of the electronic structure of matter, on 
scales ranging from single atoms to extended solids. Its relatively low computa- 
tional cost makes it attractive as a tool for answering material-specific questions 
in fields as diverse as quantum chemistry, materials science, and many-body 
physics [T|. In practice, DFT is most commonly used in its spin-dependent 
self-consistent-field formulation as Kohn-Sham spin-dcnsity-functional theory 
(SDFT) |21 El- Alternative formulations that use other variables in addition to 
(or instead of) the spin densities, or that proceed without invocation of Kohn- 
Sham-type equations, are, however, occasionally useful. Some such alternative 
formulations are described in this contribution. Section El gives a brief overview 
of three recent developments of this type in SDFT, relating to currents, non- 
collinear spin states and nonuniquencss. In each case a surprising feature of 
the respective extension of DFT will be highlighted. Section El contains a more 
detailed description of a recent application of SDFT to model Hamiltonians of 
the Heisenberg type, where results can be obtained without solving Kohn-Sham 
equations. 

2 Some recent developments in spin-density-func- 
tional theory 

From a practical point of view, among the most important tasks remaining in the 
development of SDFT are to obtain ever better approximations for the exchange- 
correlation (xc) functional, and to learn how to calculate further quantities of 
physical and chemical interest from the output quantities of SDFT, the spin- 
resolved charge densities. While total energies and related quantities are readily 
accessible in terms of these densities, many other interesting quantities are not. 

2.1 Electronic currents 

Electronic currents are one example. Such currents are functionals of the charge 
density, but cannot be calculated as such because this functional is unknown. 
The current calculated from the Kohn-Sham (KS) orbitals of SDFT is that of 
the auxiliary noninteracting system, and not guaranteed by the basic theorems 
of SDFT to be the correct current. The full physical (gauge invariant) current 
is j(r) = j p (r) + (e/mc)n(r)A(r), where A(r) is the vector potential andj p (r) is 
the so-called paramagnetic current. The paramagnetic current calculated from 
the KS orbitals of SDFT is 

jfV) = ^T^K(r)V^ fe (r)-^(r)V^(r)] ^j p [n](r). (1) 

k 

Current-density functional theory (CDFT) 01 E] provides a way to calculate the 
true equilibrium currents by obtaining the orbitals from a modified KS equation 
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that features an xc vector potential A xc (r), chosen such that in terms of these 
orbitals the above inequality becomes an equality. To achieve this, approximate 
xc functionals of CDFT must depend on the current density, in addition to the 
charge and spin densities. This dependence is a complicating factor both in 
the construction of approximate functionals and in their implementation. In 
spite of these complications, different ways of constructing CDFT functionals 
and solving the CDFT KS equations have been explored for atoms, molecules 
and solids @l El El EI El El CHI CU As an alternative to solving the full CDFT 
equations, usable when current-dependent effects are expected to be small, it 
has been proposed to cast the CDFT KS equations in the form of SDFT KS 
equations plus a remainder depending on the xc vector potential A xc (r), which 
is then treated perturbatively [T^]. 

An unexpected feature of CDFT is that, although SDFT and CDFT arc 
independent formulations of DFT for magnetic systems, there is a, somewhat 
obscure, link between CDFT and SDFT, which may be used to extract infor- 
mation on functionals of one from the other. In Ref. J3| this link is obtained 
mathematically, but its physical origin is simple to understand: the existence of 
spin currents of the form j s = V x m, where m is the spin magnetization, im- 
plies a connection between the formalisms of SDFT and CDFT, because being a 
manifestation of the spin degrees of freedom these currents must be describable 
via SDFT, while as currents they can also be built into the formal framework 
of CDFT. Some consequences of the resulting connections and consistency rela- 
tions between SDFT and CDFT are explored in [T3]. For more details on CDFT 
we refer to Refs. H El El 13 El El EB El H2 HSj 



2.2 Noncollinear magnetism 

The standard spin-up spin-down formulation of SDFT is incapable of describing 
spin configurations that are eigenstates of the square of the total spin operator 
S 2 , but not its z— component, S z . Such noncollinear configurations can arise, 
e.g., in the presence of spin-orbit coupling, domain walls, external noncollinear 
fields, or from spontaneous symmetry breaking. Noncollinear ground states 
are found, e.g., in rare-earth based magnetic systems. In magnetic excitations, 
noncollinearity is ubiquitious. Consequently, noncollinear formulations both of 
static J3JE1E1 and time-dependent ^2] DFT have been proposed, but appli- 
cations are limited to the local-spin-density approximation (LSDA). In trying 
to describe noncollinear spin states with generalized-gradient approximations 
(GGAs) or other beyond-LSDA functionals an interesting problem arises. 

Due to its origin in the electron liquid, the exchange-correlation energy in 
the LSDA depends only on the modulus of the spin magnetization m(r). The 
corresponding conjugate field, the xc magnetic field 

is then always parallel to m. Beyond the LSDA, however, the xc magnetic 
field ceases to be parallel to the spin magnetization. Mathematically this oc- 
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curs because the existence of spatial derivatives of any component of m in the 
functional leads to additional terms in the derivative that are not parallel to 
m ^1 Intuitively this can be understood because the existence of spa- 
tial gradients of the spin density introduces another preferred direction into the 
problem. In general, the xc magnetic field can thus be decomposed as 

B ;cc (r)=By c (r)+B^(r), (3) 

where the second contribution is locally perpendicular to the spin magnetization 
m(r). This perpendicular xc magnetic field exerts an xc torque oc m(r) x B xc (r) 
on the spin distribution. In the time-dependent case some consequences of this 
torque were described in JS] , but it is present even in the ground state as soon 
as noncollinearity appears and is treated beyond the LSDA. This subtle effect 
is little explored in present-day SDFT. 

Standard implementations of noncollincar SDFT assume parallelity between 
B IC and m(r), and can thus not be used consistently with GGA-type function- 
als. What the best way is to treat noncollinearity beyond LSDA is an open 
question at present. As an alternative to formulations in terms of the spin mag- 
netization, it has been proposed, e.g., to employ the spin-offdiagonal element 
of the density matrix (the so called 'staggered density' [201 ) ^ but this requires 
the construction of functionals of this new variable, and the solution of more 
complicated KS equations |21| . 



2.3 Nonuniqueness 

At first sight more remote from numerical applications are fundamental issues 
of SDFT, such as the recently discussed question of nonuniqueness of the SDFT 
potentials. The question to what extent the effective potentials of SDFT (and 
other generalized DFTs) are uniquely determined by the densities arises because 
contrary to what is sometimes claimed in the literature, the standard proofs of 
the Hohcnberg-Kohn theorem (by contradiction j22j and by constrained search 
|23| ) guarantee only the existence of a one-on-one mapping between densities 
and wave functions, not densities and potentials. The extension of the mapping 
to potentials requires, as an additional step, inversion of Schrodinger's equation 
to express the potential in terms of the universal operators T and U and the 
wave function |2S1|^]. For charge-only DFT one has simply 

P. (4) 

which shows that the ground-state wave function W (which is known to be a 
functional of the density) determines the potential V uniquely up to an additive 
constant (here the ground-state energy, E ). Since one always has an additive 
constant free in the definition of a potential, this shows that the physically 
relevant part of V, in particular its variation in space, is fully determined by Wo 
and hence by n(r). Already in the early literature on noncollincar SDFT it was 
pointed out that this inversion is problematic if one works with 2-component 



4 



spinors and has four potential components (v, B x , B y , B z ) to determine. Von 
Barth and Hedin gave an explicit example of the resulting nonuniqueness of the 
SDFT effective potentials, albeit one that only works for a single electron |25) . 

Recently, examples have been given that show that nonuniqueness is not a 
pathological feature of single-electron systems, or limited to SDFT, but a general 
phenomenon, expected to occur whenever one works with more than one density 
variable |2(:il I27| . In Ref . |27| it was shown how a large class of nonunique poten- 
tials is related to the existence of certain types of constants of motion (systematic 
nonuniqueness) while others arise from special features of the ground state (acci- 
dental nonuniqueness) . A few simple examples from collinear SDFT were given, 
illustrating the main features of both types of nonuniqueness. The simplicity of 
these examples notwithstanding, their existence shows that there cannot be a 
unique mapping between densities and potentials in generalized DFTs, in addi- 
tion to the one between densities and wave functions. More complex examples 
of nonuniqueness were found to be possible in noncollinear situations [2E1I2E], 
on lattices or in CDFT |3U]. The implications of nonuniqueness for practi- 
cal SDFT calculations are still under study H3 HO HO Ell 132] One conclusion 
that has emerged from this body of work is that the additional mapping between 
density and potential must be regarded as a special feature of spin unpolarized 
(charge-only) DFT. Surprisingly, the HK theorem of unpolarized (charge-only) 
DFT is thus considerably stronger than its counterparts in, e.g., SDFT and 
CDFT. 

3 Spin density-functional theory for the Heisen- 
berg model 

Open problems, such as the best way to deal with currents, noncollinearity 
and nonuniqueness, provide fascinating challenges, but they have never stoped 
DFT from advancing on the practical side. In a typical application of DFT one 
deals with the ab initio many-electron Hamiltonian comprising kinetic, poten- 
tial, and (Coulomb) interaction terms. By means of the Hohenberg-Kohn (HK) 
theorem and the Kohn-Sham construction this many-body problem is mapped 
on an auxiliary noninteracting one. The search for good approximations for 
the exchange-correlation potential appearing in the auxiliary KS Hamiltonian 
is one of the most active fields of DFT. However, the HK and KS techniques 
are not restricted to the ab initio Hamiltonian. We can simply regard them as 
efficient tools for mapping some many-body problem on a simple one-body prob- 
lem. In many areas of science, most notably perhaps in statistical physics and 
many-body physics, but occasionally also in quantum chemistry, the fundamen- 
tal Hamiltonian is replaced by a simpler model Hamiltonian before embarking 
on any numerical calculation. The list of such simplified Hamiltonians is long, 
and comprises important models such as those of Ising, Heisenberg and Hub- 
bard. This section is devoted to a description of a recent reformulation and 
application of SDFT to model Hamiltonians of the Heisenberg type [551 EH] , 
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3.1 The Heisenberg model 



Introduced in 1928 by Dirac and Heisenberg an d discussed in some detail 

as a model for ferromagnetic solids in 1932 by van Vleck |37], the Heisenberg 
model is the paradigmatic model for the magnetism of local magnetic moments 
Sj. In quantum chemistry it has, e.g., been used in the description of the 
magnetic properties and spin states of conjugated hydrocarbons [33 EH and a 
large variety of metallic and organometallic compounds 0U1 BS 1421 14"5| . 
In its simplest form the Heisenberg model reads 



also have an orbital contribution. The sums run over nearest neighbours among 
the N sites of a lattice of, in principle, arbitrary geometry, and J parametrizes 
the nearest-neighbor interactions. We see qualitatively that a parallel (ferro- 
magnetic) or antiparallcl (antifcrromagnetic) alignment of Sj and Sj is favored 
according to whether J is negative or positive. 

In the presence of an external magnetic field Bj the model takes the form 



If Bj and Sj = (Sj) are the same everywhere, the model is homogeneous, i.e., all 
its sites are equivalent. In all other cases it is inhomogeneous. Note that even 
for Bj = (no external field) inhomogeneous states are possible. Such solutions 
may arise because the self-consistent many-body ground state need not preserve 
all symmetries of the Hamiltonian, or because the Hamiltonian itself comprises 
inequivalent lattice sites. The latter is the case, e.g., in systems with impurities, 
or, more generally, with of different types of magnetic atoms. 

In spite of its formal simplicity, exact analytical results for the Heisenberg 
model are known only for homogeneous linear chains with S = 1/2, by means of 
the Bethe ansatz ^J05]. For higher spin, higher dimension, or inhomogeneous 
situations, exact solutions are known only numerically, for small systems (less 
than « 40 sites), where the model can be diagonalized exactly by Lanczos 
techniques. In the presence of inhomogeneities, such as one substitutional spin 
of value S' ^ S, no analytical solution is known. 

The simplest approximation to the Hamiltonian (jSJ) is the mean-field (Hartree- 
like) approximation, in which the vector operators S are substituted by classical 
vector spins S. The mean- field ground-state energy per spin, in units of J, is 
thus easily obtained for a homogeneous infinite system of spins S in d dimensions 
(with periodic boundary conditions): 




(5) 



where the S, are angular- momentum vector operators satisfying S 2 \Sm) = 
h 2 S(S + l)\Sm). In accordance with common terminology we refer to the ex- 
pectation values of Sj, the magnetic moments Sj, as spins, although they may 




(6) 



e™ F (S,d) 




(S,d) 



-dS 2 , 



(7) 



NJ 
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where d is the dimensionality of the linear (d — 1), square (d = 2) or cubic 
(d = 3) lattice. Other lattice geometries can be dealt with in the same way. An 
improvement on the mean-field estimate for eo(S, d) is obtained from spin- wave 
theory |46| . according to which 

e sw {S: d) = E^S,d) = _ dg2 + d _ 1/5 ^ _ ^ s (8) 

Here we used the scaling hypothesis conjectured in Ref. 133], which accounts 
for the numerical spin- wave results to three decimal places, to write the re- 
sults in different dimensionalities in closed form, as a function of d. Beyond 
spin- wave theory, sophisticated numerical many-body techniques, such as the 
density-matrix renormalization group (DMRG) 07J or Quantum Monte Carlo 
simulations |48| . have been employed to further improve the ground-state en- 
ergy. 

The resulting expressions for E$, however, suffer from a major limitation, 
namely, they are applicable only to homogeneous systems, in which all lattice 
sites are equivalent. This limitation is restrictive, because in real problems it 
is very common to find inhomogeneities, induced for instance by external non- 
uniform magnetic fields, anisotropic crystal-fields, finite-size effects, boundary 
conditions, local impurities, lattice defects, etc. These rather common situa- 
tions introduce very significant difficulties in traditional approaches like Bethe 
ansatz, DMRG, spin-wave theory or even Monte Carlo, by breaking transla- 
tional symmetry. In order to make progress in dealing with inhomogeneous 
Heisenberg models we need a many-body technique capable of handling spatial 
inhomogeneity in large systems. Density-functional theory is our first choice in 
this regard. 

3.2 Spin-distribution functionals and local-spin approxi- 
mation 

Density-functional theory |2] offers, with the local-density approximation 
(LDA), a simple prescription for obtaining ground-state properties of an inho- 
mogeneous system, based on the knowledge of a related homogeneous system. 
While in ab initio applications of DFT this homogeneous system is the electron 
liquid, in our context it is the homogeneous Heisenberg model, in which all sites 
are equivalent. 

Ref. [33J shows that the basic ingredients of DFT (in particular the Hohenberg- 
Kohn theorem 22 ) stilll hold for Hamiltonians of the form @, provided one 
uses instead of the ground-state density n(r) the ground-state expectation value 
of the spin vectors, Si = (\E , |Sj|\f r ). The Heisenberg-model Hohenberg-Kohn the- 
orem then states that the ground-state expectation value of any observable O 
is a functional of the distribution Si, i.e., 0[Sj] = (*[Sj]|0|*[Sj]). In this 
Heisenberg-model DFT the above mean-field expressions take the place of the 
Hartree term in Coulomb DFT, and the correlation energy e c (S) is defined as 
the difference eo(5 l ) — e^f F (S), where 60(5*) is the full ground-state energy. From 
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the above, the spin- wave approximation for e c (S) is thus 

e sw = e sw _ 6 mf = d -i/ 5 (~-l)S. (9) 

The Heisenberg- model Hohenberg-Kohn theorem, proved in |33| , guarantees 
that in inhomogeneous systems the total energy is a functional of the spatial 
distribution of the classical vectors S^. Hence 

E = EolSt] = £ Mi? [S,] + £ C [S,] . (10) 

Since the mean- field term is already a functional of Si, we just need an approxi- 
mation for the functional i? c [S*]. This is obtained using a local-spin approxima- 
tion (LSA), patterned after the conventional LDA of Coulomb DFT. The LSA 
consists in replacing locally the variable S by Si in e c (S): 

N 

E c [S t ] « E e [Si] LSA = £ e e (S)\s^\ Bt \ ■ (11) 

i 

Using the above spin-wave approximation for e c (S) we obtain, e.g., 

E L3A-SW [Shdj jj = d -l/5j \Si\ , (12) 

^ n ' i=i 

and thus the ground-state energy functional in this approximation is 

More sophisticated expressions, based on fully numerical evaluation of 60(6*) 
are also available |.''>.'ij . but in first applications the simple approximation \Y6\ 
turned out to be already reasonably accurate |34| . 

As always in dealing with LDA-type functionals one must distinguish two 
sources of error. One is the LDA itself, the other the quality of the descrip- 
tion of the underlying homogeneous reference system. In Coulomb DFT several 
alternative parametrizations of the LDA are available |251 , which differ con- 
siderably in form, and slightly in results. In the present context, too, one can 
construct various LSA functionals, depending on the level of sophistication em- 
ployed in solving the uniform Heisenberg model. For example, for an infinite 
linear chain of spin S — 1/2 we have from the above spin- wave based LSA 
e /{JN) = -0.25 - 0.181690 = -0.431690, which is 2.6% from the exact Bethe 
ansatz value —0.443147. Interestingly, we find this same margin of deviation also 
when the LSA is used for inhomogeneous situations (see below). This seems to 
indicate that the LSA concept as such is quite reliable for Heisenberg models. 

From the above, it should have become clear that the formalism of DFT for 
the Heisenberg model is built in complete parallelity to that of Coulomb DFT 
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(or that for the Hubbard model [SOI)- One interesting difference is that within 
spin-wave theory the dimension dependence of the functional is approximately 
known. This is very different from the ab initio case, in which nothing is known 
about this dependence. In general, however, one may expect that LSA for the 
Heisenberg model shares the virtues and defects of its close cousin, the LDA for 
the Coulomb problem: efficient access to large and inhomogeneous systems on 
one side, offset by moderate accuracy on the other. In the next subsection we 
demonstrate these aspects by analysing two interesting inhomogeneous models 
within LSA. (Other applications of the above functionals can be found in Ref. 

E1-) 

3.3 Applications to inhomogeneous Heisenberg models: 
boundaries and impurities 

The usual DFT prescription for obtaining ground-state energies is to solve Kohn- 
Sham equations. These equations, however, were originally [2] introduced to 
deal with the kinetic-energy term in the Coulomb Hamiltonian. The Heisen- 
berg model does not have such a term, and direct minimization of the en- 
ergy functional is much simpler than indirect minimization via self-consistency 
(Kohn-Sham) equations. In practice we just need to minimize the total energy 
of the system, written as a functional of the spin distribution, with respect to 
Si. This already represents a considerable improvement over the mean- field ap- 
proximation with no extra computational cost; minimization of Eq SA ~ sw [Si] 
is no more complicated than that of Eq IF [Si] . 

Lets first consider a finite ring of N = 16 spins S = 1/2, with periodic 
boundary conditions. The mean field approximation yields eo = —0.25, while 
the LSA — SW approximation gives eo = —0.43169, about 3% from the exact 
value —0.44639 obtained by direct diagonalization. Let us now replace any spin 
of the ring by an impurity spin Si. This slight modification is enought to chal- 
lenge any analytical approach, but since the ring is not too big, it is still possible 
to numerically diagonalize H in order to get the exact ground state. Figure U 
shows how the ground-state energy depends on Si, and compares the values ob- 
tained by mean-field and LSA-SW calculations with the numerically exact ones. 
Up to Si = 2.5 the deviation from the exact values are of order 3%. Above 
this value the deviation increases, showing that the LSA approximation is not 
good for large impurity spin (corresponding to a very rapid variation of local 
magnetic properties). For impurity spins up to Si = 5/2, on the other hand, 
LSA-SW is seen to provide a substantial improvement on the mean-field ap- 
proximation. Unlike the exact calculations, LSA-SW can easily be extended to 
systems with more impurities and more sites. It may thus provide a convenient 
way of estimating (and improving on) the error of the mean-field approximation 
in more complicated systems. 

Next we consider 4x4 square lattices with an impurity of spin Si at the cor- 
ner, at one side, or in the interior of the lattice. Here we use open boundary con- 
dition and thus even without impurity the systems are inhomogeneous. Results 
are collected in Tab. ^ The first column is for Si = 1/2, i.e., for the uniform 
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system, and the value listed differs 6.6% from the exact value eo = —0.574325 
obtained from numerical diagonalization; we believe the values for other impu- 
rity spins are within this accuracy, although no exact values are available, as 
far as we know. Besides the obvious energetic tendency of higher impurity-spin 
to prefer the interior of the lattice to its surface, which is already visible in the 
mean-field results, these data also incorporate less trivial correlation effects: In 
mean- field theory, for example, a spin 3/2 at the corner of the square is degen- 
erate with a spin 1 in the interior of the lattice (e^f F = —0.4375 in both cases). 
This degeneracy is broken by correlation effects. Similarly, mean-field theory 
predicts degeneracy between a spin 5/2 at one side and a spin 2 in the interior 
( e <f F = —0.5625 for both cases). Again this degeneracy is removed by including 
correlations. 

4 Summary 

Kohn-Sham spin-density functional theory is doubtlessly the most widely used 
formulation of density-functional theory. However, in spite of its great popu- 
larity and usefulness some interesting fundamental questions still remain open. 
The existence of spin currents, providing a link between SDFT and CDFT; 
the existence of an exchange-correlation torque appearing in noncollinear spin 
configurations; and the nonuniqueness of the potentials of SDFT and other gen- 
eralized DFTs are examples of the kind of surprise the SDFT formalism still 
holds. 

Such open questions notwithstanding, the utility of (S)DFT is extending 
even into areas remote from ab initio calculations. Recognition that the basic 
tools and concepts of DFT — such as the Hohcnberg-Kohn theorem, and the 
local-density approximation — are not restricted to the original Coulomb prob- 
lem allows one to apply these tools and concepts also to many other interesting 
inhomogeneous model Hamiltonians. In particular, density-functional theory 
together with results from mean-field and spin-wave theory provide a simple 
manner to obtain estimates of the ground-state energy for spin distributions 
described by the inhomogeneous Hciscnberg model. This is here illustrated by 
results for a finite ring and a 4 x 4 lattice, both with substitutional impurity 
spins Si — 1,3/2,2 or 5/2. While ring calculations may find applications in 
modelling the spin states of hydrocarbons, the square-lattice data illustrate a 
way in which LSA can be useful in nanomagnetism: To predict the structure of 
self-assembled magnetic nanostructures it is clearly important to know which 
lattice sites are degenerate in mean-field theory, and which of these degeneracies 
are removed by correlations. 

This work was supported by FAPESP and CNPq. We thank P. E. G. Assis, 
and F. C. Alcaraz for useful discussions and collaboration on the Heisenberg 
model. K.C. thanks L. N. Oliveira, E. K. U. Gross and G. Vignale for use- 
ful discussions and collaborations on various aspects of spin- density- functional 
theory. 
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Tabic 1: LSA-SW ground-state energy —E /NJ of a 4 x 4 square lattice of 
5=1/2 sites with a substitutional impurity of spin Si at the corner (four 
equivalent positions), on a side (eight equivalent positions) or in the interior of 
the lattice (four equivalent positions). 
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Figure 1: Per-site ground-state energy of a ring of N = 16 spin S =1/2 sites as 
a function of the impurity spin Si. The LSA results, based on spin-wave theory, 
(dashed line) deviate from the exact ones (circles) by about 3%. The mean field 
approximation (continuous line) deviates by up to 50%. The inset illustrates 
the geometry under study, for the case of N = 10. 



